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ABSTRACT 

We present high spatial resolution Submillimeter Array observations and supplemen- 
tary single-dish photometry of the molecular gas and dust around IRAS 04158+2805, 
a young source with spectral type M5-M6 in the Taurus star-forming region. A bright, 
highly elongated dust structure that extends 8" (~1120 AU) in diameter is revealed in 
a 883 /im thermal continuum image. This emission geometry is in good agreement with 
optical observations that show a similar structure in absorption, aligned perpendicu- 
lar to bipolar scattered light nebulae. However, the interferometric data also clearly 
demonstrate that the submillimeter continuum emission is not centrally concentrated, 
but rather appears to have a toroidal geometry with substantially lower intensities inside 
a radius of ~250-300AU. Spatially resolved emission from the CO J=3— 2 transition 
exhibits a velocity gradient along the major axis of the dust structure. If this kinematic 
pattern is interpreted as the signature of rotation around a central object, a relatively 
low mass is inferred (M* ~ 0.3 Mq, with a ~50% uncertainty). We discuss several possi- 
ble explanations for the observed gas and dust environment around IRAS 04158+2805, 
including a flattened envelope with an outflow cavity and a large circumbinary ring. 
This source offers unique views of the gas and dust environment surrounding a young 
low-mass stellar system. Its properties are generally not commensurate with formation 
scenarios for such low-mass objects that rely on dynamical ejection, but rather con- 
firms that a single mechanism — molecular cloud core collapse and fragmentation — can 
produce stars over a wide range of stellar masses (at least an order of magnitude). 

Subject headings: circumstellar matter — stars: formation — stars: low-mass, brown 
dwarfs — stars: pre-main-sequence — stars: individual (IRAS 04158+2805) 
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Introduction 



Circumstellar material profoundly influences the star formation process. A large-scale envelope 
acts as the local mass reservoir during the collapse and growth of a central protostar; a more compact 
disk regulates how that material is transported onto the star itself. A great deal of progress has 
been made in understanding the physical conditions present in the gas and dust surrounding Sun- 
like stars (M* ~ 0.5-3 M©). However, it is unclear if those results can be extrapolated across the 
entire stellar mass spectrum. The low end of this spectrum (M* < 0.3 Mq) is of particular interest, 
as observations of the material around young low-mass objects^ can help resolve an ongoing debate 
over their dominant formation mechanism. 

The observed abundance of such low-mass objects relative to their more massive counterparts 
is difficult to account for with the traditional model for isolated star formation. Two distinct 
modifications have been proposed: (a) the tu rbulent fragmentation of cloud cores into smaller 
building units (e.g., iPadoan fe Nordlundl 12004 ) : and (6) a dynamical alternative where one com- 
ponent of a multiple system is ejected and p rematurely cut off from its accretion reservoir (e.g., 
Reipurth &: Clarke! l200ll : lUmbreit et al.1 120051 ) . Both scenarios make distinct and conflicting pre- 
dictions about the gas and dust environments that should be associated with low-mass objects. 
In the former, a disk/envelope structure similar to those noted around higher mass T Tauri stars 
would be expected. And for the latter, only a t runcated disk (w ith an initial Rd < 10 AU) would 
survive the tidal stripping and ejection process ([Bate et al.ll2002l ). I f the remnant material is suffi- 
ciently viscous, it could spread to radii up to ~100AU in ~lMyr (jArmitage Clarkdll997l ). but 
would have shed a vast majority of its initial mass in the ejection process. Therefore, observational 
constraints on the structure of the gas and dust around low-mass objects offer one avenue to help 
distinguish how they form. 

A variety of observations demonstrate that the signatures of circumstellar gas and dust at 
small ra dii (up to a few AU ) are indeed common for young low- mass objects (see the recent re- 
view by iLuhman et al.l 120071 ) . More detailed individual studies confirm that such material has a 
geometry and composition similar to the disks around more m assive T Tauri stars (jNatta Testi 



2001 



Pascucci et al 



(Klein et al 



2003 



Scholz et al 



200< 



2003 :lAllers et al.ll2006l : iBouy et al.ll2008l ). although typically at lower masses 



However, these unresolved observations can not unambigu- 
ously constrain some key properties of this material, most significantly its spatial structure and 
extent. 

In this article, we present new submillimeter observations of the molecular gas and dust 
surrounding IRAS 04158+2805. Located in the ~1 Myr-old Taurus star-for ming region, IRAS 
0415 8 +2805 has a cool central so urce with a n estimated sp e ctral type of M5-M6 (jWhite fc Hillenbrand 



2004 



Luhman 



2006 



Beckl 120071 ). Recently, iGlauser et al.1 (|2008l ) provided an initial analysis of this 



For simplicity, we will refer to objects with M, < 0.3 M as low-mass objects, including both very low-mass stars 
and brown dwarfs. 



- 3 - 



source and the dust structure surrounding it by modeling the broadband spectrum along with op- 
tical/infrared scattered light images. They suggested that their observations are best explained 
by a large (diameter of 2240 AU), massive (dust mass of l-2xl(r 4 M ; gas/dust ratio of 220±J^) 
circumstellar disk around a low-mass star (M* ~ 0.1-0.2 Mq), with no need to include any contri- 
bution from a more extended envelope. Our new data represent the first high spatial resolution 
view of the gas and dust environment around this and similar cool, young objects at submillimeter 
wavelengths. They also provide an opportunity for a rare, albeit crude, estimate of M* from the 
spatio-kinematics of circumstellar gas that is independent of pre-main-sequence evolution models. 
These observations are addressed in §2, and the resulting data products are highlighted in §3. In 
§4, we discuss the structure of the gas and dust surrounding IRAS 04158+2805 and what it can 
tell us about the central source and the formation of low-mass objects in general. The results are 
summarized in 85. 



2. Observations and Data Reduction 

Submillimeter interferometric observations of IRAS 04158+2805 (a = 4 h 18 m 58 f !14, 5 = +28 t T2'23"6 
[ J2000] ) were conducted in both the extended and compact configurations of the Submillimeter 



Array (SMA; iHo et al.ll2004l ) on 2006 December 8 and 2007 January 27, respectively. In these con- 
figurations, the 8 SMA antennas (6 m diameter each) spanned baselines from ~15-225m. Double 
sideband receivers with a total bandwidth of 4 GHz were tuned to an effective continuum frequency 
of 339.9 GHz (883 ^m). The correlator was simultaneously configured to observe the CO J=3— 2 
transition at 345.796 GHz at a spectral resolution of 0.70 km s -1 . 

The observations cycled between IRAS 04158+2805 and two gain calibrators (3C111 and 3C84), 
with 20 minutes on target and 10 minutes on calibrator. The data were obtained in excellent 
observing conditions, with zenith optical depths at 225 GHz of 0.05-0.07. The raw visibilities were 
calibrated with the IDL-based MIR package Passband calibration was conducted with bright 
quasars (3C273, 3C279), and complex gain calibration was performed with 3C111. Using 3C84 as a 
consistency check on the gain solution, we estimate that phase noise generates an effective "seeing" 
of <0 / /l. Titan and Callisto were used to set the absolute flux scale, which is accurate to ~10%. 
The MIRIAD package was utilized for the standard tasks of Fourier inversion, deconvolution, and 
imaging of the calibrated visibilities. Continuum and line maps from the combined datasets were 
made with natural weighting, providing a synthesized beam FWHM of 1'.'12 x / /86 at a position 
angle of 92°. 

IRAS 04158+2805 was also mapped at 450 and 850 /im using the jiggle-mode of the SCUBA 
camera at the 15 m James Clerk Maxwell Telescope (JCMT) on 2003 November 12 in dry, stable 
conditions (r ~ 0.05 at 225 GHz). Flux calibration accurate to ~10% and 25% at 850 and 450 fim, 



http : //cfa-www .harvard. edu/$\sim$cqi/mircook .html . 
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respectively, was achieved with observations of Neptune, and pointing was referenced to DG Tau 
and CRL 618. The bright emission, 0.44 ± 0.04 Jy at 850 fim and 1.6 ± 0.4 Jy at 450 ^m, was not 
resolved at either wavelength in 15" and 9" beams, respectively. An additional map at 350 fim 
was obtained with the SHARC-II camera at the 10 m Caltech Submillimeter Observatory (CSO) 
on 2004 October 3 in similarly dry weather. Those data confirm an unresolved source (9" beam) 
with a total flux density of 2.0 ± 0.5 Jy. All errors on flux densities include absolute calibration 
uncertainties (10% at 850 /im, and ~25% at shorter wavelengths). 

Optical images of IRAS 04158+2805 in the i?-band were obtained with the University of Hawaii 
2.2 m telescope on 2004 October 6. After the standard reduction of CCD data, the 7 x 300 s images 
were co-added and convolved with a 3-pixel (0 / .'7) FWHM Gaussian kernel to match the typical 
seeing for the observations and improve sensitivity to faint emission. The observations were not 
photometric. 



Results 



Images of the gas and dust surrounding IRAS 04158+2805 are shown in Figure HJ The optical 
i?-band (~0.7/mi) image in Figure QJi reveals a high aspect ratio dark lane of absorbing material 
oriented roughly E- W (shown as g r aysca le-filled contours) that can be traced nearly 8" in diameter. 
First discovered bv lGlauser et al.l (|2008l . see also Watson et al. 2007), this silhouette lane obscures 
the faint nebulosity associated with the nearby Herbig Ae star V892 Tau. Also pointed out by 
Glauser et al.l (|2008l ). the image in Figure QJi shows the conical nebulae perpendicular to the dark 
lane that are often produced by starlight scattered off the surface of flattened dust structures. 

Figure [lb shows the first high resolution image of the thermal continuum emission from IRAS 
04158+2805, at a wavelength of 883 /im. The morphology of that emission is remarkably similar 
to the optical silhouette, extending roughly 8" across, corr esponding t o a diameter o f 1120 AU at 
the distance of the Taurus clouds (d as 140 pc: lEliaslll978l ). Note that I Glauser et al.l (|2008l ) trace 
/-band scattered light for this source to twice this diameter (~16", or 2240 AU). The difference 
compared to the i?-band image in Figured^ is related to the comparatively poor sensitivity of our 
observations, as was also noted for the infrared scattered light images in their study. Because the 
submillimeter emission was clearly unresolved in the 350 and 450 /im maps, we can infer that it 
has a diameter of < 9" (1260 AU). The absence of very short antenna spacings implies that the 
interf erometric observations have significantly diminished sensitivity to structures larger than ~10" 
(e.g., IWilner &: Welch! Il994l ). Therefore, because of the diminished sensitivity to faint submillime- 
ter emission on large scales, there is no reason to suspec t that the sized infer red here is physically 
inconsistent with the scattered light image presented by iGlauser et al.l (|2008l ) . The integrated flux 
density at 883 is 407 ± 41 mJy, consistent with the SCUBA measurement. A fit of the SMA 
visibilities with an elliptical Gaussian model indicates an inclination of 62 ± 3° (where 90° corre- 
sponds to edge-on) at a position angle of 93 ± 1° (measured east of north), in ex cellent agreement 
with the constraints imposed by modeling the scattered light (iGlauser et al.ll2008l ). Fits with a ring 
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or pedastal emission model give the same results. 

The high surface brightness at large distances from the central source and large spatial extent 
of the continuum emission are fairly unique. The steep drop of the visibilities over a short range of 
baseline lengths, as shown in Figured! indicates that essentially all of the emission is concentrated 
at large spatial scales (radii > 250 AU). More significantly, the visibilities show a pronounced null 
at ~30-40kA and very little flux on longer baselines. These are definitive signatures of a brightness 
distribution with a central "hole" , inside of which there is significantly diminished intensity. Assum- 
ing a simple ring g eometry, the location of the null suggests an inner radius of ~250-300 AU (e.g., 



Hughes et al.l 120071 ). Asymmetries in the submillimeter map confirm this geometry, with emission 
peaks centered roughly ±1'.'9 (~265 AU) from the image center. The western lobe is approximately 
50% brighter than the eastern lobe. These asymmetries are real, as they are clearly detected in 
3 separate observing runs at the same locati on and intensity (the 2 deta iled here, and the hybrid 
compact / extended observations described by Andrews Williams 2007a ] E 



A velocity- integrated intensity (0 th moment) map of the CO J=3— 2 emission from IRAS 
04158+2805 is shown in Figure [It. The CO emission morphology is significantly different than the 
optical silhouette and submillimeter continuum, occupying a region ~4" on a side with a central, 
slightly elongated peak. The integrated intensity is 40.7+0.4 Jy km s -1 (using 9 channels for a total 
width of 6.3 km s _1 ). Fainter CO emission extends above and below the plane of the continuum 
emission in an "X" -shaped pattern out to ~2" from the image center. There is a steep drop in 
the CO intensity along the plane of the silhouette that corresponds precisely to the positions of 
the submillimeter continuum peaks described above. The kinematic structure of the CO associated 
with IRAS 04158+2805 is exhibited in Figure [3l These channel maps clearly show a CO velocity 
gradient along the major axis of the continuum emission and silhouette (i.e., in the east-west 
direction), similar to the pattern expected for rotation around a central source. 



4. Discussion 



4.1. The Central Object 



Estimating the masses of young stars and brown dwarfs generally requires reference to theoret- 
ical models of their structural evolution. The quantitative reliability of such models for individual 
stars at ages <1 Myr is highly uncertain, particularly for the low end of the stellar mass spec- 
trum where complicated convection physics, dusty atmospheres, and uncer tain molecular opacities 



and initial conditions can strongly affect the obse rvational diagnostics (e.g.. iD'Antona &: Mazzitelli 



1994 



Baraffe et al.l 120021 ; iMontalban et al.l 12004) . Dynamic al constraints on M* from the orbital 



properties of either a companion star (e.g., iMathieul Il994l ) or gas in a circumstellar disk (e.g., 



3 In fact, the image made from the hybrid configuration in lAndrews fc Williams 
peak of this asymmetry, and heavily resolves the remainder of the emission. 



(|2007al ) shows only the western 
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Simon et al.ll2000l ) provide an extremely valuable indepen de nt check that c an po tentially be used 



to calibrate these models (see Hillenbrand White! 12004 ). IStassun et al.l (|2006l ) report the only 
such measurements of young ~M6 dwarfs to date, for an eclipsing binary system in Orion. With 
the first spatially resolved measurements of molecular gas in apparent rotation around a young low- 
mass object with similar spectral type, the CO observations of IRAS 04158+2805 presented here 
can provide another such estimate of M#, although admittedly with a great deal more uncertainty. 

Figured] shows the position-velocity diagram of the CO J=3— 2 line emission for IRAS 04158+2805, 
where the angular offset is the distance from the observed phase center (such that positions to the 
east have positive values) and the velocity offset is taken relative to the systemic value (V* = 7.4 km 
s _1 ). Overlaid on the diagram is an inclined (i = 62°) Keplerian rotation profile that matches 
the general kinematic trend in the data rather well (thick curve), with a central point mass of 
M* = 0.30 Mq. Additional rotation profiles for M* = 0.15 and 0.45 M Q are also shown for refer- 
ence. It should be noted that these profiles assume pure Keplerian orbital motion around a central 
source, which may not be the case in the more complex IRAS 04158+2805 environment. For now, 
we assume that the dominant kinematic trend is such rotation, but potential complications will be 
addressed below. 

A comparison of this dynamical M* estimate with predictions from stellar evolution models is 
a challenge because of uncertainties in the luminosity (due to scattered light contamination at short 
wavelengths), extinction, and effective temperature for the central source. Recent analysis of the 
scattered light spectrum f rom IRAS 04158+2805 suggest that this source is cool, with initial spectral 



type estimates of M6+1 (White fc Hillenbrand 



2004 



Beck] [20071) refined to M5.25+0.25 (optical) 



and M6+0.5 finfrared: EuhmaJ20od ')FI Using the latter classifications and the lLuhman et al.l (|2003l ) 
empirical effective temperature scale, we estimate T* = 3050 + 75 K. Unfortunately, large uncertain- 
ties i n the extinction, where estimates range from Ay ~ 9-16 (e.g. JWhite fe Hillenbrandll2004l ; iBeck 
20071 ). make a luminosity determination difficult. However, we can use pre-main-sequence models 



for an assumed age in the inferred tempera ture range to estimate a centra l mas s. Using the range 
quoted above and an age of ~lMyr, the b'Antona fc Mazzitellil (|l997h and baraffe et al.l dl998l ) 



models indicate M* ~ 0.09-0.16 M Q . This mass range does not include the additional uncertainties 
on the adopted effective temperature scale or the unknown age of the source. 

The M* values inferred from the pre-main-sequence models lie roughly a factor of 2 below 
the nominal value that best describes the kinematics of the CO gas inferred above. The unquoted 
and poorly understood uncertainties in both measurements may be the simple explanation for this 
apparent discrepancy. Without a detailed physical model for the local gas structure around IRAS 
04158+2805, it is difficult to estimate a statistically acceptable range of M* from the CO kinematics; 
there may very well be overlap with the higher end of the range from the pre-main-sequence 
models. The effective temperature scale adopted above is also uncertain. A ~500K increase in the 
effective temperature (to a value that is typically associated with M2 stars) would reconcile the 



4 Note that previous spectral type assignments ranged from ~K7 to M3 ( Kenvon et al. 



19981 ; 



Luhman 



19991 ). 
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M» es timates. Indeed, using an independent method of fitting infrared spectra, iDoppmann et al 
l|2005l ) infer T* = 3500 K for the IRAS 4158+2805 central sou r ce, le ading to M* ra 0.35-0.45 M 
in the aforementioned models. Moreover, iHillenbrand &: White! ([2004]) have shown that pre-main- 
sequence models generally tend to underpredict M* compared to dynamical measurements in this 
low-mass range, although typicall y only by < 2 %. In their recent study that explicitly includes 
the effects of scattered starlight, iGlauser et al.l (|2008l ) argue for a very luminous central source, 
w 0.4 Lq. While this value must be somewhat degenerate with the dust properties and structure 
assumed in their model, it is clearly the most robust estimate available. The pre-main-sequence 
models imply that such a luminous cool source would be very young, with an age significantly less 
than 1 Myr. The model mass tracks at such ages are not well understood in light of the uncertain 
initial conditions. 

Clearly, the aforementioned uncertainties can be solely responsible for the apparently differ- 
ent M* estimates from the CO kinematics and pre-main-sequence models. Alternatively, there 
is also an appealing and simple physical explanation — the IRAS 04158+2805 central source 
may be a roughly equal-mass binary. In this scenario, the M* estimate from the CO kinemat- 
ics refers to the total mass of the stellar system (i.e., the mass interior to the gas that produces 
the line emission), while the estimate from the pre-main-sequence models a ssumes only a single 
star. Therefore we would expect to see both a high luminosity (as claimed bv lGlauser et al.l l2008) 
and M*(CO) ~ 2JVf* (m odels) . A variety of surveys suggest that s uch binaries are fairly common, 
~35- 45% for M dwarfs [Fischer M arcv!ll992l:l Reid Gizislll997h and ~10-30% for cooler objects 
(e.g., iBurgasser et al.1 120071 ). While it is difficult to definitively state the properties of the central 
source(s), we can further assess the nature of IRAS 04158+2805 by analyzing the structure of the 
material that surrounds it. 



4.2. The Gas and Dust Environment 



Given the elongated structures shown in Figure [IJ it is reasonable to suggest th at the gas and 



dust a round IRAS 04158+2805 reside in an exceptionally large circumstellar disk. IGlauser et al 



(|2008l ) have recently demonstrated that they are able to fit high-quality scattered light images and 
the SED for this source with a simple disk structure, and without the need to invoke any envelope 
component. With the new data presented here, it is worthwhile to re-examine the structure of this 
material. 

IRAS 04158+2805 is exceptionally bright at submillimeter wavelengths, with an 850um lu- 
minosity larger than ~85% of all other sources detected in the lAndrews &: Williams! (|2005l ) Tau- 
rus survey. Breaking that sample down further, the IRAS 04158+2805 submillimeter emission is 
brighter than ~90% of Class II sources (disk only), but is comparable to the median Class I source 
(disk + envelope). If we adopt the standard emissivity for disks (0.034 cm 2 g" 1 at 883 /mi, includ- 
ing a gas-to-dust ratio of 100) and assume that the su bmillimeter continuum emission is optically 



thin and has a characteristic temperature of ~20 K (e.g.. lBeckwith et al. 



1990 



Andrews &z Williams 
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20051 . l2007bl ) , the integrated SMA flux density suggests a large total mass of gas and dust is present 
around IR AS 04158+280 5 , ^0. 03 M^. This value is in good agreement with the disk mass range 
inferred bv lGlauser et al.l ([20081 ). ~ 0.02-0.04 Mq. Estimating masses for circumstellar material 



in this way is inherently uncertain (perhaps by an order of magnitude), due to the challenge of 
observationally constraining the optical properties of dust grains and the assumed mass conversion 
of a trace constituent (dust) to the presumably dominant species (molecular gas). 

Perhaps more compelling (and easier to interpret) is the unusual geometry that is observed 
for the dust s tructure around IRA S 04158+2805. The enormous size of the scattered light neb- 



ulae noted by iGlauser et al.l (|2008l . 2240 AU diameter) and the submillimeter continu um emission 



presented here (1120 AU diameter) is exceptionally rare for circumstella r disks (e.g., iPietu et al 



2007 



Watson et aI1l2007l : kndrews fc Williams! l2007al: Irlughes et al.lboosl). rivaled only b v the 



cir- 



cumbinary disk around UY Aur (jDuvert et al.lll998l : IClose et al.lll998l : IPotter et al.ll2000l ). On the 

contrary, the observed size is fairly typical (if not on the small end) for larger-scale envelopes (e.g., 



Looney &: Mundy 



2000 



Eisner et al. 



200.51 ). 



Moreover, our resolved observations of the submillimeter continuum reveal that a large central 
region (perhaps ~500 AU in diameter) has significantly diminished intensity at 883 /im, suggesting 
a ring-like or toroidal geometry for the dust structure. Similar geometries have been inferred for 
a handful of circumstellar disks, and a variety of underlying causes are possible. The so-called 
"transition" disks have their inner regions (out to a few tens of AU in radius) largely cleared of 
observable material, presumably by disk evolution processes li ke particle growth, photoevaporation , 



or dynamical inter actions with a young planetary system (e.g. JPietu et al. 



2006 



Hughes et al. 



2007 



Brown et al.ll2008l ). However, given the large spatial scale of the central depression and the absence 
of a corresponding dip in the infrared part of the IRAS 04158+2805 SED, as shown in Figure 
this scenario is unlikely. 

Alternatively, extremely dense disks with edge-on orientations can potentially render even the 
subm illimeter emissio n optically thick, resulting in a central depression in its brightness distribution 
(e.g., IWolf et alJl2008l ). Although the IRAS 04158+2805 dust structure is not oriented edge-on, 
perhaps longer pathlengths through the exceptionally large structure can provide compensating 
column densities. If the bulk of the submillimet er emission is optically thi ck, the spectrum should 
have a spectral index a = 2, where F v oc u a (jBeckwith fc Sargentl Il99ll ). Power-law fits to the 
submillimeter photometry are significantly steeper than the optically thick case, with a = 2.7 + 0.3 
(0.45-1.3 mm; a = 3.8 ± 0.6 for 0.85-1.3 mm). At these wavelengths, such a steep spectrum is only 
noted for a few percent of Class II sourc es in Taurus and p Oph, but is m ore common (25-35%) 



in the envelopes around Class I objects (lAndrews Williams 



20051 . l2007bl ). The observed steep 



submillimeter SED would be diffic ult to produce unl ess a large fraction of the dust was optically thin 
at those wavelengths. As noted by I Wolf et al.l (|2008l ). a clear test of this possibility can be made with 
resolved continuum observations at A > 883 //m. In such data, the material should become optically 
thick only at significantly smaller radii than at 883 ^m, resulting in a more centrally concentrated 
emission profile (or, equivalently, the null in the visibilities should be detected at a comparatively 
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longer deprojected baseline). Note that if the material inside this central depression is optically 
thick at 883 /xm, an enormous reservoir of mass is not accounted for in the above estimate. 

A third possibility for explaining the submillimeter continuum depression is that the emitting 
dust resides in a circumbinary ring. The material around close binary stars is expected to be 
dynamically cleared on scales similar to the physical separation of the stellar components (e.g., 
Jensen et al"1ll996l : iGuilloteau et al.1 ll999h Fl Numerical simulations of this dynamical interaction 
suggest that a circumbinary disk gap would be cleared out to a radius ~1.5-3x the projected 
semimajor separation of the stellar components, and the individual circu mstellar disks would be 
truncated at a radius ~0.2-0.5x the projected semimajor separation (e.g.. lArtvmowicz Lubow 



1994 ) . If the 883 /zm continuum depression noted here is caused by such clearing, we would expect 
the projected binary separation to be ~90-180 AU (~0 / /6-l / /3) for a reasonable range of eccentricities, 
and the remnant individual circumstellar disks might be truncated at radii of ~20-90AU. This 
configuration is similar to the UY Aur system, where a ring with ~20" diameter is detected in both 
scattered light and CO around a roughly equal- mass b inary with a similar projected separation 



(Duvert et al 



1998 



Close et al.l 119981 : iPotter et al.1 1200G ) . The primary difference is in the masses 
involved; the stellar mass in the IRAS 04158+2805 system is ~4x less than for UY Aur, but the 
circumstellar mass (or rather 850 /im luminosity) is ~4x higher. While this possibility is certainly 
appealing in its reconciliation of both the information on the central source (see §4.1) and the 
structure of the gas and dust emission, testing the multiplicity of IRAS 04158+2805 will be a 
challenge given the obscuration of the central region by the dust structure. Decomposing high 
resolution spectra of scattered starlight may represent the best opportunity in this case. 

The IRAS 04158+2805 SED, shown in Figure exhibits a strong infrared excess with a flat or 
slightly rising slope from 2-25 /im. A variety of solid-state absorption features are noted in higher 
resoluti on infrared spectra, including H2O ice (3 /im; not shown), silicates (10 /mi), and CO2 ice 
(15 /mi; iBeckl 120071 ; iFurlan et al.ll2008l ). The bright and steep submillimeter spectrum has already 
been discussed above. While all of these SED properties are not necessarily inc o nsistent with a large, 
cold, highly inclined circumste llar disk (e.g., iMen'shchikov &; Henningjll997l : IChiang &: Goldreich 
1999J; iPontoppidan et al.ll2005l ). taken together they are more common for C l ass I/Flat-Spectrum 
objects that harbor a remnant accretion enve lo pe (e.g.. IWhitnev et al.1 120031 : IWatson et al.1 [2004]; 



Boogert et al.1 120041 : IPontoppidan et al J 120071 ). IFurlan et al.1 (120081 ) demonstrate that the IRAS 



04158+2805 SED features are strikingly similar to those for several other Class I sources in Taurus 
(see their Fig. 2) and can be successfully reproduced in detail with an envelope model (although 
the new observations of this source warrant a modification of their geometric assumptions). So, all 
of the observational evidence presented here is also consistent with a flattened envelope structure 
around IRAS 04158+2805. Such a flattened envelope geometry is predicted by molecular cloud 



^Recent observations suggest that the circumbinary scenario may be the more appr opriate explanation in some 



cases for the diminish ed central emission in both transition (jlreland fc Kraus 
( Guilloteau et alJbOQgh . 



2008f ) and dense, edge-on disks 
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core collapse models that incorporate rotation or magnetic fields (jTerebev et al.lll984l : iGalli fe Shu 
19931 ). and have been clearly detected in other cases, albeit on significantly larger spatial scales 
than noted here (e.g.. lLoonev et al.ll2007l ). 



In this scenario, a partially evacuated outflow cavity at the center of t he flattened envelope 
may explain the central depression in the sub millimeter continuum data (e.g.. IWhitney et alJl2003i . 



their Models 3 or 4). Idauser et all (|2008l ) indicate that IRAS 04158+2805 drives a substan- 



tial Ha jet oriented in the north-south direction, perpendicular to the observed dust structure. 
The velocit y-integrated CO J=3 — 2 morphology does resemble the bases of some molecular out- 



flows (e.g., lArce fc Sargentl 120061 ) . but the absence of both a velocity gra dient in the proposed 



flow direction (nort h-south) and larger scale molecular outflow signatures (jBontemps et alJll99a : 
Gomez et al.lll997l ) are difficult to reconcile with that interpretation. Nevertheless, some contri- 
bution to the CO emission from an outflow and even the envelope (i.e., motions unrelated to 
Keplerian rotation) would complicate the observed spatio-kinematics. The dynamical estimate of 
M* discussed above should be treated with caution with regards to this possibility. Of course, the 
envelope and binary ideas are not mutually exclusive, and dynamical clearing of the inner envelope 
and disks is still a viable alternative. 

The images in Figures [Tfc and [3] demonstrate that, unlike the submillimeter continuum, the 
CO emission is centrally concentrated; most of the line emission is from radii inside the toroidal 
structure traced in the continuum. While this may suggest that the bulk of the line emission 
is generated in the disk(s) interior to that structure, it is difficult to provide a reliable origin 
without a more complete physical model for the IRAS 04158+2805 environment. Future studies 
of this source should focus on searching for definitive molecular outflow signatures, reso lving dense 



gas tracers that are commonly associated with envelopes (e.g., lJ0rgensen et al.l 120071 ). observing 
CO transitions at higher spectral resolution to provide a better tracer of the velocity field, and 
reconciling the optical/near- infrared scattered light observations with the submillimeter results in 
an effort to converge on a consistent model of the local gas and dust structure. 

Considering the low stellar mass range implied for IRAS 04158+2805, it is difficult to reconcile 
the large mass reservoir and size of this environment with an y formation scenario that relies on it 
being dynamically ejected from a more massive system (e.g 



Reipurth Clarke 



2001 



Bate et al 



20021 ). Rather, the observed gas and dust structure is more consistent with the standard scenario, 
where an accretion disk and perhaps a remnant envelope developed around a single (or binary) 
source during the gravitational collapse of a molecular cloud core fragment. Regardless of the 
precise value of M*, IRAS 04158+2805 is a remarkable example of the fact that the standard 
picture for star formation worked out for higher-mass stars is also applicable at the low end of the 
stellar mass function. 
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5. Summary 

We have presented new submillimeter observations of the molecular gas and dust around IRAS 
04158+2805, a young low-mass object in the Taurus star-forming region. The data reveal a complex 
and surprising environment, from which we infer the following: 

1. A high aspect ratio (i « 62°) dust structure with a projected diameter of ~1120 AU and 
estimated (gas+dust) mass of ~0.03 M is noted in optical absorption and submillimeter continuum 
emission images. The resolved 883 fim data clearly demonstrate a central depression in the emission, 
suggesting a ring-like geometry with an inner radius of ~250-300 AU. A variety of structures that 
may explain the observations are investigated, with a flattened envelope or circumbinary ring 
perhaps the most likely. 

2. Spatially resolved CO J=3— 2 line emission shows a velocity gradient along the major axis 
of the dust structure. Assuming the gas traced by this transition is in Keplerian orbital motion 
around a central source, the spatio-kinematic properties of the emission are used to estimate a 
central mass, M* rs 0.3 M q (with ~50% uncertainty). This estimate is roughly a factor of 2 higher 
than the mass range inferred from a reasonable effective temperature range and pre-main-sequence 
models, ~0. 09-0. 16 M , but the significant uncertainties on both measurements do not formally 
rule out their mutual consistency. However, a roughly equal-mass binary system represents an 
intriguing possibility in this case, as it would naturally account for this possible M* discrepancy, a 
very high luminosity, and the ring-like geometry of the dust structure. 

3. The large spatial extent and mass of the complex gas and dust environment surrounding 
IRAS 04158+2805 would be difficult to retain in the dynamical ejection models used to explain 
the formation of such low-mass objects. Instead, this case appears to illustrate that the standard 
gravitational collapse models applied to stars with masses up to ~2-3 M can be extrapolated down 
to objects with M* ~ 0.1-0.3 M . More detailed models of the IRAS 04158+2805 gas and dust 
structure and specialized observations are recommended. 

We thank F. Menard and G. Duchene for encouraging our interest in this source, and Elise 
Furlan for kindly providing the mid-infrared spectrum. We are grateful to an anonymous referee 
for a helpful, constructive review and for emphasizing the potential circumbinary nature of IRAS 
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Institution and the Academia Sinica. Support for this work was provided by NASA through 
Hubble Fellowship grant #HF-01203.01-A awarded by the Space Telescope Science Institute, which 
is operated by the Association of Universities for Research in Astronomy, Inc., for NASA, under 
contract NAS 5-26555. MCL and KNA acknowledge support for this work from NSF grant AST- 
0407441. 
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Fig. 1. — Images of the gas and dust surrounding IRAS 04158+2805. (left) Optical i?-band (0.7 /xm) 
image showing a dark lane in silhouette (grayscale contours; linear steps) against background 
nebulosity and scattered starlight (unfilled contours; logarithmic steps), (middle) Submillimeter 
(883 /xm) dust continuum image from the SMA. Contours start at 3 a (~3.5mJy beam -1 ) and 
increase by factors of \/2. (right) Velocity-integrated (0 th moment) emission image from the CO 
J=3— 2 transition, constructed from 9 spectral channels (a velocity width of 6.3 km s -1 ). Contours 
start at 3 a (0.5 Jy km s -1 beam -1 ) and increase by factors of \pl. Effective FWHM (synthesized) 
beam sizes are shown in the lower right corner of each panel. 
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Fig. 2. — The real part of the IRAS 04158+2805 883 /jm visibilities as a function of the interferom- 
eter baseline, a fter deprojection and elliptical averaging according to the observed geometry (see 
Lav et alJll997l ). The steep drop in the visibilities at short baselines and the pronounced null at 
~30-40kA confirm the morphology in Fig. lb; the submillimeter emission is not centrally peaked. 



-17- 




4 2 -2-4 
Ac, ["] 



Fig. 3. — Channel maps of the IRAS 04158+2805 CO J=3-2 line emission. Contours start at 3 a 
(0.2 Jy beam -1 ) and increase linearly in 6 a increments. The FWHM synthesized beam sizes are 
shown in the lower right of each panel, and the LSR velocity of each channel (in km s -1 ) is shown in 
the upper left. The crosshairs mark the horizontal extent of the submillimeter continuum emission 
and the vertical extent corresponding to the inclination (62°) derived from fits of the continuum 
visibilities. There is a clear velocity gradient along the major axis of the continuum emission and 
optical silhouette. 
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Fig. 4. — Position-velocity diagram of the CO J=3— 2 line emission. Contours start at O.Uy 
and increase in steps of \[2. The ordinate shows the projected offset from the image center (with 
positive values to the east), and the abscissa shows the velocity offset from the systemic value, 
+7.4 km s _1 . The heavy curve represents an inclined (i = 62°) Keplerian rotation profile that is 
most similar to the data, for = 0.3 M Q . Lighter curves on either side show similar profiles for 
M* = 0.15 and 0.45 M . 



-18- 



BO 
O 




10 100 
\ [jum] 



1000 



Fig. 5. — The IRAS 041 58+2805 spectral energy distribution. Optical photome try was taken from 



Luhman et al 



Furlan et al 



Luhman fc Riekd (119981 ) . in frared data from the 2MASS Point Source Catalog (jCutri et al.ll2003l ) 



(2006|), and iBeichman et al.l (|1986l ). the mid-infrared Spitzer IRS spectrum from 



(200 8J), the submillimeter data from measurements presented here, and the 1.3 mm 



flux density from lMotte & Andrei (fcoOll ). See iFurlan et al.l (|2008l ) for a more detailed discussion of 
the mid- infrared spectral features. 



